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Abstract 

We analyse the nonlinear Kuramoto-Sivashinsky equation to de- 
velop an accurate finite difference approximation to its dynamics. The 
analysis is based upon centre manifold theory so we are assured that 
the finite difference model accurately models the dynamics and may 
be constructed systematically. The theory is applied after dividing the 
physical domain into small elements by introducing insulating inter- 
nal boundaries which are later removed. The Kuramoto-Sivashinsky 
equation is used as an example to show how holistic finite differences 
may be applied to fourth order, nonlinear, spatio-temporal dynamical 
systems. This novel centre manifold approach is holistic in the sense 
that it treats the dynamical equations as a whole, not just as the sum 
of separate terms. 



1 Introduction 

We continue exploring the new approach to finite difference approximation 
introduced by Roberts jllj by approximating the dynamics of solutions to the 
Kuramoto-Sivashinsky equation |3], ^ [5|. In some non-dimensional form we 
take the following partial differential equation (pde) to govern the evolution 

of u(x, t): 

du du j^ 2u ^ u q n\ 
dt dx dx 2 dx A 
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1 Introduction 



This model equation includes the mechanisms of linear growth u xx controlled 
by the parameter R, high-order dissipation, u xxxx , and nonlinear advec- 
tion/steepening, uu x . Consider implementing the method of lines by dis- 
cretising in x and integrating in time as a set of ordinary differential equa- 
tions. A finite difference approximation to (0) on a regular grid in x is 
straightforward, say Xj = jh for some grid spacing h. For example, the 
linear term 

d 2 u Uj + i — 2uj + Uj-i 



dx 2 h 2 

However, there are differing valid alternatives for the nonlinear term uu x : 
two possibilities are 
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Which is better? The answer depends upon how the discretisation of the 
nonlinearity interacts with the dynamics of other terms. The conventional 
approach of considering the discretisation of each term separately does not 
tell us. Instead, in order to find the best discretisation we consider the 
influence of all terms in the equation in a holistic approach. 



As introduced in (10| and discussed in §0, centre manifold theory [Q. ||, 
e.g.] has appropriate characteristics to do this. It addresses the evolution of 
a dynamical system in a neighbourhood of a marginally stable fixed point; 
based upon the linear dynamics the theory guarantees that an accurate low- 
dimensional description of the nonlinear dynamics may be deduced. For 
example the analysis herein supports the first approximation in (fj) but with 
higher-order and nonlinear corrections. The analysis of the Kuramoto-Siva- 
shinsky equation (|I|) in § [3] favours the discretisation 

Uj(-u j+2 + 9u j+ i - 9ii,-_i + Uj-z) 



duj 
— - + 
dt 



+ 



16h 

'U 2 +1 - u]_ x \ (U j+2 U j+1 - Uj-. 2 Uj-l 



+ R 



16h J V 48/i 

-Uj+2 + 16Mj + i — 30Uj + lQUj-i — Uj^2 

12^2 



Uj+2 - 4u j+1 + 6Uj - 4Mj_i + Uj-2 _ n ,o\ 

+ h4 -U, [6) 

as a low-order approximation. Provided the initial conditions are not too 
extreme, centre manifold theory assures us that such a discretisation models 
the dynamics of (|]) to errors O (||m|| 3 , h 2 ). Such accuracy on a relatively 
coarse grid is extremely useful for such stiff pdes. Further, because the 
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centre manifold is composed of actual solutions to the dynamical system, we 
are assured that equation @ models the whole of the Kuramoto-Sivashinsky 
equation, independent of its algebraic form. 

The discretisation (H) is a low-order approximation, centre manifold the- 
ory also provides systematic corrections. Analysis to higher orders in the 
nonlinearity, discussed in § [3], shows higher order corrections to the discreti- 
sation of the nonlinear terms. The specific finite difference models presented 
here were derived by a computer algebra program. Computer algebra is an 
effective tool because of the systematic nature of centre manifold theory ||. 

In this preliminary exploration of the approximation of the Kuramoto- 
Sivashinsky equation (P we only consider an infinite domain or strictly pe- 
riodic solutions in finite domains. Then all elements of the discretisation 
are identical by symmetry and the analysis of all elements is simultaneous. 
However, if physical boundaries to the domain of the pde are present, then 
those elements near a physical boundary will need special treatment. Further 
research is needed on this and other issues. 



2 Centre manifold theory underpins the fi- 
delity 

Here we describe in detail one way to place the discretisation of the Kura- 
moto-Sivashinsky equation (|1|) within the purview of centre manifold theory. 

The discretisation is established via an equi-spaced grid of collocation 
points, Xj = jh say, for some small spacing h. Here we scale the Kuramoto- 
Sivashinsky to the scale of the grid spacing h. Thus our view of the dynamics 
shrinks with h. This is different to the analysis in ]IU| and allows the linear 
dynamics to be dominated by simply the highest order spatial derivative 
term. We work on the scale of the grid by transforming ([!]) to the following 
space and time scales: £ = x/h and r = t/h A , giving 



du , o du , 9 „<9 2 w d 4 u 
— + fru— + h 2 R-— + — - 
Or d£ de <9£ 4 



0. 



(4) 



Then the crucial step: at midpoints £j+i/2 = (£j+£j + i)/2 artificial boundaries 
are introduced: 



du + du 

d 3 u+ _ d 3 u~ 









ir 

av 
at 2 



u 

d 2 u 



(5) 
(6) 
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2 Centre manifold theory underpins the fidelity 



where u + is just to the right of a midpoint and u to the left. The introduc- 
tion of the near identity operator 

^ 1 + 5-i + ^U... = ^coth(^, (7) 
12 720 30240 2 \ 2 / 

ensures that high-order approximations to linear terms are obtained exactly 
as discussed in |L0], §4]: it is remarkable that the exactly equivalent operator 
works for both Burgers' equation and the Kuramoto-Sivashinsky equation. 
These boundaries divide the domain into a set of elements, the jth element 
centred upon £j and of width A£ = 1. A non-zero value of the parameter 
7 couples these elements together so that when 7 = 1 the grid scaled Ku- 
ramoto-Sivashinsky pde (^) is effectively restored over the whole domain by 
ensuring sufficient continuity between adjacent elements. 

The application of centre manifold theory is based upon a linear picture 
of the dynamics. Adjoin the dynamically trivial equations 

and consider the dynamics in the extended state space («(£), 7, /i). This 
is a standard trick used to unfold bifurcations (l|, §1.5] or to justify long- 
wave approximations ||. Within each element u = 7 = is a fixed point. 
Linearized about each fixed point, that is to an error O {\\u\\ 2 + 7 2 + h 2 ), the 
pde is 

du d 4 u du 
s.t. 



dr <9£ 4 ' <9£ 



£=±1/2 



d 3 u 



= 

£=±1/2 



namely the hyperdiffusion equation with essentially insulating boundary con- 
ditions. There are thus linear eigenmodes associated with each element: 

U(X [ eXnT cos[mr(£ - &-1/2)] , < f < £ j+1/2 , (g) 

' 1 , otherwise , 

for n — 0, 1, . . ., where the decay rate of each mode is A n = — n 4 7r 4 ; together 
with the trivial modes 7 = const, h = const and u = 0. In a domain with 
m elements, evidentally all eigenvalues are negative, — 7r 4 or less, except for 
m + 2 zero eigenvalues: 1 associated with each of the m elements and 2 from 
the trivial (||). Thus, provided the nonlinear terms in (|j) are sufficiently well 
behaved, the existence theorem (p|, p281] or |Tl], p96]) guarantees that a 
m + 2 dimensional centre manifold M. exists for The centre manifold 

M. is parameterized by 7, h and a measure of u in each element, say uf 
using u to denote the collection of such parameters, M. is written as 

u(Z,T) = v(Z;u, 1 ,h). (10) 
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In this the analysis has a very similar appearance to that of finite elements. 
The theorem also asserts that on the centre manifold the parameters Uj evolve 
deterministically 

Uj = 9j(u,l,h) , (11) 

where iij denotes duj/dr, and gj is the restriction of (||-^|) to M.. In this 
approach the parameters of the description of the centre manifold may be 
anything that sensibly measures the size of u in each element — we simply 
choose the value of u at the grid points, Uj{j) = u(£j, r). This provides the 
necessary amplitude conditions, namely that Uj = v(^j; u,j, h). 

The above application of the theorem establishes that in principle we may 
find the dynamics (|TTD of the interacting elements of the discretisation. A 
low order approximation written in unsealed variables is given in (|3|). 

The next outstanding question to answer is: how can we be sure that such 
a description of the interacting elements does actually model the dynamics of 
the original system Here, the relevance theorem of centre manifolds, 

|2|, p282] or Jll], pl28], guarantees that all solutions of which remain in 



the neighbourhood of the origin in 7> h) space are exponentially quickly 

attracted to a solution of the m finite difference equations ([Tl]). For practical 
purposes the rate of attraction is estimated by the leading negative eigen- 
value, here — 7r 4 . Centre manifold theory also guarantees that the stability 
near the origin is the same in both the model and the original. Thus the finite 
difference model will be stable if the original dynamics are stable. After ex- 
ponentially quick transients have died out, the finite difference equation (|TTD 
on the centre manifold accurately models the complete system 

The last piece of theoretical support tells us how to approximate the shape 
of the centre manifold and the evolution thereon. Approximation theorems 
such as that by Carr & Muncaster |2|, p283] assure us that upon substituting 
the ansatz (|10|-|TT|) into the original and solving to some order of error 

in 7 and h, then A4 and the evolution thereon will be approximated to 
the same order. The catch with this application is that we need to evaluate 
the approximations at 7 = 1 because it is only then that the artificial internal 
boundaries are removed. In some applications of such an artificial homotopy 
good convergence in the parameter 7 [|7| has been found. Thus although the 
order of error estimates do provide assurance, the actual error due to the 
evaluation at 7 = 1 should be also assessed otherwise. Here we have crafted 
the interaction ([]) between elements so that low order terms in 7 recover the 
exact finite difference formula for linear terms. Note that although centre 
manifold theory "guarantees" useful properties near the origin in («(£), 7, h) 
space, because of the need to evaluate asymptotic expressions at 7 = 1, we 
have used a weaker term elsewhere, namely "assures". 
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3 Numerical comparisons show the effectiveness 



3 Numerical comparisons show the effective- 
ness 

We now turn to a detailed description of the centre manifold model for the 
Kuramoto-Sivashinsky equation ([!]). 

The algebraic details of the derivation of the centre manifold model (|TU|- 
TTD are handled by computer algebra.Q In an algorithm introduced in || , the 



program iterates to drive to zero the residuals of the governing differential 
equation (f|) and its boundary conditions (0-0). Hence by the Approximation 
theorems we assuradly construct appropriate approximations to the centre 
manifold model. 

The finite difference model is given by the evolution on the centre mani- 
fold. In order to represent the spatial fourth derivative in the Kuramoto-Siva- 
shinsky equation we need to determine the interactions between next-nearest 
neighbouring elements. Thus the first approximation we can consider involves 
quadratic terms in 7. After returning to x and t variables it is 

duj 7.R , n . 7 . . 

— = -— (u j+ i-2uj+Uj-i)- — Uj(u j+1 - Uj-i) 

7 2 
7 2 i? 

+ (u j+2 u j+1 + 3u j+2 Uj - 3u j+1 - 3u j+1 Uj 

+3Uj-\Uj + 3Uj_i — 3lLj~2Uj — Uj^2Uj-l) 

l 2 h 2 r.a 



H Ua (— 30ui+2 — 170-u 1+1 + 256m,- — 170uji — 30ujo) 

60480 L 3 

+ Uj (—12Quj + 2Uj + i — 54uj + iUj_i — 12GiLj^2 u j-i) 
+ u^ +1 (10uj +2 — 20uj + i + 235^) 
+w 2 „ 1 (10^_i - 20u i _i + 235uj) 

+ 0(|| M || 4 , 7 3 ,/ i 4 ) • (12) 

The first two lines recorded here, when evaluated for 7 = 1, form the conven- 
tional second-order finite difference equation for the Kuramoto-Sivashinsky 
equation ([!]). The third line when evaluated for 7 = 1 gives the fourth order 

1 The reduce computer algebra source code is available from the authors upon request. 
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Figure 1: Contours of an accurate solution u(x, t), , to compare with 



numerical approximations (|T^) : • • •, the conventional approximation, errors 
O (h 2 ); — — — , the first correction, errors O (||w|| 3 ); — • — • — , the second 
correction, errors (||ii|| 4 ). Kuramoto-Sivashinsky equation fl]]) with param- 
eter R = 2 is discretised on just m = 8 elements in [0, 2ir) and drawn with 
contour interval Am = 3. 



accurate corrections to the Ru xx term. The fourth and further lines above 
start accounting systematically for the variations in the field u within each 
element and how they affect the evolution through the nonlinear term. Fi- 
nite difference equations derived via this approach holistically model all the 
interacting dynamics of the entire pde. 

To show the effectiveness of the approach we compare the finite difference 
model obtained from various truncations of ( |P2"D to accurate numerical solu- 
tions obtained on a much finer grid. Choosing m intervals on [0, 2tc) gives 
an element length h = 27r/m and grid points Xj = jh for j = 0, . . . , m — 1. 
Because of the antisymmetry in u(x, t) about x = kir, when starting from the 
initial condition u(x, 0) = 10 sin x, we only display the interval [0, 7r]. There 
are three different approximations from (0), with 7 = 1, depending upon 
where the expansion is truncated. The first two lines form a model with 
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O (h 2 ) errors (a conventional finite difference approximation), the first five 
lines provide the first correction with (9(||u|| 3 ) errors, and all shown terms 
form the model with 0(||w|| 4 ) errors. The solutions of these models over 
< t < 1 with m = 8 and R = 2 are shown in Figure [l|. Observe that the 
leading approximation (dotted) is significantly in error whereas the next two 
refinements (dot-dashed and dashed) are overall more accurate, especially 
near the peak. Such accuracy is remarkable considering the nonlinearity, 
and the few points in the discretisation, m — 8. 



4 Conclusion 

Centre manifold theory is a powerful new approach to deriving finite differ- 
ence models of dynamical systems. Many details need to researched for a 
general application of the theory. However, there are many promising fea- 
tures of this application to the Kuramoto-Sivashinsky equation (|l|) and the 



earlier example of Burgers' equation [ HI . 
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